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Heavy quark potential and quarkonia dissociation rates 
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Abstract. Quenched lattice data for the QQ interaction (in terms of heavy quark free energies) in the 
color singlet channel at finite temperatures are fitted and used within the nonrelativistic Schrodinger 
equation formalism to obtain binding energies and scattering phase shifts for the lowest eigenstates in 
the charmonium and bottomonium systems in a hot gluon plasma. The partial dissociation rate due to 
the Bhanot-Peskin process is calculated using different assumptions for the gluon distribution function, 
including free massless gluons, massive gluons, and massive damped gluons. It is demonstrated that a 
temperature dependent gluon mass has an essential influence on the heavy quarkonia dissociation, but 
that this process alone is insufficient to describe the heavy quarkonia dissociation rates. 



PACS. 12.38.Gc Lattice QCD calculations 
S=C=B=0, mass > 2.5 GeV 



12.38.Mh Quark-gluon plasma - 14.40.Gx Mesons with 



1 Introduction 

Heavy quarkonia have been suggested as hard probes of 
the quark gluon plasma since the modification of static 
interactions at finite temperature eventually implies a dis- 
solution of heavy quarkonia bound states into the contin- 
. £h ' uum °f scattering states (Mott effect). This effect results 
, in a suppression of heavy quarkonia production in heavy- 
ion collisions as an observable signal. Since the Mott tem- 
' peratures for J/tp, T and T 1 as obtained by solving the 
Schrodinger equation for a screened Cornell-type poten- 
tial lie well above the critical temperature T c for decon- 
fincment 2] , it has soon been realized that a kinetic theory 
is necessary for the description of heavy quarkonia dissoci- 
ation PJ, see [4"||^] for recent formulations. Solutions of the 
Schrodinger equation provide the basis for the evaluation 
of cross sections and rates for the Bhanot-Peskin process 
of heavy quarkonia dissociation by gluon impact . In 
this contribution we present a new fit to the singlet free 
energies from quenched lattice QCD simulations 8 , 9 and 
show that binding energies and cross sections deviate from 
those obtained for Debye potential fits [TUIllllll2| . 



2 Heavy quark potential 

The main source of our knowledge of the quark-antiquark 
static interaction at high temperature originates from cal- 
culations of the free energy for a quark-antiquark system 



from the Polyakov-loop correlator from lattice QCD. The 
color averaged correlator receives contributions from the 
color singlet and color octet channels 
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(Tr[i(0)]Tr[it( r )]) = iex P (-^M 

(1) 

The color-singlet part F\ is extracted from the equation 

(Tr[i(0)Lt(r)])=exp('-^i) (2) 

and the color octet part is found by subtraction. The free 
energy of a quark-antiquark system thus obtained is inter- 
preted hereafter as the effective static interaction potential 
for a quark-antiquark pair surrounded by gluons. 



2.1 Zero temperature 

As there exist no color octet mesons in the vacuum, we 
will use the color singlet free energies as a potential in the 
Schrodinger equation description of the heavy quarkonia 
spectrum. Two additional parameters are to be defined 
in fitting the spectra: the heavy quark masses and a con- 
stant shift of the whole potential. The behavior of the 
quark-antiquark interaction in the color singlet channel 
was investigated in Ref. within the combined lattice 
and perturbative QCD approach. We fit these data points 
implementing a \ 2 minimization to the Ansatz 
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Fig. 1. Lattice data for the singlet free energy 9 at different 
temperatures, see legend. Solid lines are obtained with the fit 
on the basis of the Dixit formula, dashed line is the T = 
potential. 



where F% > s h ort (r) describes the interaction at short dis- 
tances whereas i ? i j i ong (r) is responsible for the long-dis- 
tance forces. Both expressions are matched at r — tq that 
is defined below. This point, as we shall see, lies in the 
domain of perturbative QCD. We use the combined linear 
and Coulomb potential to describe the long-distance inter- 
action and the Coulomb interaction with the r-dependent 
coupling constant a(r) for short distances 



4o(r) 
3 r 
16tt 

a(r) = —— 



33 v>(r 2 /c 2 ) 



r 2 — c 2 
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The formula for a(r) is obtained by solving the one- loop 
renormalization group equation for the running coupling 
constant in QCD followed by the pole subtraction |14j . 
The constant c^fa w 1.816 and the point r^^fa w 0.031 are 
determined in units of the string tension ^fa = 0.42 GeV 
from the condition that the potential is a smooth function 
at r = r*o . The result is given by the dashed line in Fig.^ 



2.2 Singlet free energy at high temperature 

For the singlet free energy of a static quark-antiquark sys- 
tem at high temperature we use the Ansatz 



Fi(r,T) = 



^l.shortM, r<r 0l 
Fi.\on g (r,T), r>r , 



(7) 



where for the short-range interaction governed by pQCD 
the zero-temperature form of the previous subsection is 
assumed. 

The long-distance interaction F± t i ons (r, T) requires the- 
oretical assumptions about its shape. Instead of the fre- 
quently used screened Coulomb potential, we follow Dixit 
[T5*| and assume the potential behavior at large r as V ~ 
exp(— ^ir 2 )/^/r. Above T c , the expression for the potential 
without the Coulomb interaction was deduced in ^B] 
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and contains two parameters: q and fx. A similar param- 
eterization was used by Digal et al. ^B]. We add to this 
term a Coulomb attraction at short distances and obtain 



Fi,ion S (r, T > T c ) = F^ ong {r, T > T c ) + F ljShoit e 

(9) 

The parameters q and /i are determined from a fit to lattice 
data the results are shown in Fig.Q] 



3 Quarkonia in a hot gluon plasma 

3.1 Quarkonia at zero temperature 



The masses of quarkonia in the vacuum are defined as 



M = 2m Q + E + v 



(10) 



where tuq is the quark mass, and the energy E is an eigen- 
value of the Schrodinger equation 



-V 2 /™ Q + V(r)]i/>(r) = E^{r) , 



(11) 



where the potential V(r) is identified with the zero tem- 
perature free energy Fi(r) of Eq. J3J up to an unknown 
constant vq. Substituting the wave function iprdm(f, 0, (f>) = 
r _1 R n i(r)Yi m (9 ', <fi) into (|X 1|) - one obtains an equation for 
R n i(r). At large r, the potential is linear, and the solution 
of this equation behaves as the Airy function R n i(r) ~ 
Ai(ftr — £), where k 3 — tjiqu and ^ = itiqE/k 2 . 

The masses of IS and 4S states are used as input. For 
charmonium we obtain m c = 1.45 GeV and the constant 
vo = —302 MeV. For the bottomonium we have m\, = 
4.785 GeV with the same vq. Once m c , mb, and vq are 
fixed, the remaining quarkonia spectrum is described |17j . 



3.2 Quarkonia at finite temperature 

The Schrodinger equation for a bound state in the QGP 
has the form 

[-V 2 + m Q V cS (r,T)]il>(r,T) = m Q E B (T)i;(r,T),(12) 

where Eb(T) > is the temperature dependent binding 
energy. The medium effects on the QQ system are mod- 
eled using the singlet free energies as an effective potential 
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Fig. 2. Binding energies of heavy quarkonia states from solu- 
tion of the Schrodinger equation for the T-dependent effective 
potential of Fig. Q 



V cS (r, T) = Fx(r, T) - F^T) with the continuum thresh- 
old Foo(T) = linv^oo Fi(r, T). The temperature depen- 
dent mass of a quarkonium bound state is defined as 
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M(T) = 2m Q - E B (T) + v + F^T) 



(13) 
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Fig. 3. Scattering phase shifts for heavy charmonia and bot- 
tomonia states from solution of the Schrodinger equation for 
the T-dependent potential of Fig. 



The solutions for the binding energy both for charmonium 
and bottomonium are shown in Fig. [21 

The wave function of an unbound quark-antiquark sys- 
tem can be calculated via the S-wave phase shift function 
3s( r ) by solving the equation |18| 



dSs(k,r,T) _ m Q V cS (r,T) 



dr 



k 



sm(kr + 5 s (k,r,T))} 



(14) 

The phase shift is defined as Ss(k,T) = Ss(k,oo,T) and 
results are shown in Fig. [3] for charmonia and bottomonia 
states at different temperatures. In accordance with the 
Levinson theorem, the scattering phase shift at threshold 
changes by tt once a bound state merges the continuum 
at its Mott temperature; T Mott /T c = 1.05, 1.20, 2.25 for 
T', J/-0, T, respectively, see Figs. 2 and 3. 



4 Dissociation of quarkonia by gluon impact 

We calculate the cross section for the quarkonium disso- 
ciation after a gluon impact similarly to the calculation of 
the deuteron decay via the photon absorption |11U19| 



AnagQ (fc 2 + kp 
a (QQ)g\ UJ ) = — o 1 I / uip{r)u is {r)r dr 



(15) 

\u is (r)\ 2 dr = 1, a gQ = a s /6, k 2 = m Q E B (T), (16) 



where we used R n i(r) — u n i(r)e k ° r . For the IP state, 
one can use the wave function of a free quark-antiquark 



system: 
u\ P (r) 



sin fcr 
kr 



cosfcr, k 2 — ttiqIu) — Eb{T)). (17) 



For the constant a s in (|16fl . we take an average over the 
low energy region below 1 GeV, which gives a s « 0.48. As 
a result, we obtain the cross sections shown in Fig.0] Their 
peak values correspond to the geometrical ones (ttR 2 (T)) 
with the T-dependent radius R(T) of the quarkonia wave 
function [T7] . 

Now we can estimate the dissociation rate of the char- 
monium and bottomonium by gluon impact according to 



r (QQ) 9 ( T ) = I dsA ( s ) 
o 



duj 2 /— 
4^^ 



sa (QQ)g( UJ ) n a( UJ )^ 
(18) 



where A(s) is the normalized gluon spectral function 



A(s) 



V«7 



7r (s — m 2 ) 2 + S7 2 ' 



dsA(s) 



(19) 



and the thermal gluon distribution function is given by 
n g (uj) = 2(N 2 - l)[exp(w/T) - l] -1 . The temperature 
dependent gluon mass m g and damping width 7 are taken 
from a recent fit to lattice QCD data for the entropy in 
pure gauge [2U| . 



m 2 = 2na T 2 , 7 



3<5Tln(2.67/a), 



(20) 



where a = 12tt/[33 ln(3.7(T - T s )/T c ) 2 ] with T s = 0.46 
GeV. The results are shown in Fig. [5] and compared with 
the cases when the damping width and also the gluon mass 
are neglected. 
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Fig. 4. Cross sections 
Bhanot-Peskin process. 
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ing energies, entailing large dissociation cross sections for 
the Bhanot-Peskin process. The corresponding dissocia- 
tion rates have been evaluated using different assumptions 
for the gluon distribution function, including free mass- 
less gluons, massive gluons, and massive damped gluons. 
We have demonstrated that that a temperature depen- 
dent gluon mass has an essential influence on the heavy 
quarkonia dissociation. However, the Bhanot-Peskin pro- 
cess alone is insufficient to describe the quarkonium dis- 
sociation process 017]. On the basis of the spectrum and 
wave functions obtained here we will study next the QQ 
spectral functions above the Mott temperature and com- 
pare the results with corresponding lattice studies |21) . 
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Fig. 5. Dissociation rates of heavy quarkonia in a hot gluon 
plasma obtained with massive damped, massive undamped and 
massless (free) gluon distributions in Eq. 11181 . respectively. 



5 Conclusions 

We have used a new fit to recent quenched lattice data 
for the QQ singlet free energies at finite temperatures to 
obtain binding energies and scattering phase shifts for the 
lowest eigenstates in the charmonium and bottomonium 
systems within the Schrodinger equation formalism. In 
contrast to results on the basis of a Debye potential fit, we 
obtain much smaller finite temperature quarkonia bind- 
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